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ABSTRACT 

in 

\ We discuss that energetic electrons are generated near the X-type magnetic recon- 

I nection region due to a surfing acceleration mechanism. In a thin plasma sheet, the 

polarization electric fields pointing towards the neutral sheet are induced around the 
boundary between the lobe and plasma sheet in association with the Hall electric cur- 

■ rent. By using a particle-in-cell simulation, we demonstrate that the polarization electric 
Cn| : fields are strongly enhanced in an externally driven reconnection system, and some elec- 
trons can be trapped by the electrostatic potential well of the polarization field. During 

\ the trapping phase, the electrons can gain their energies from the convection/inductive 

\ reconnection electric fields. We discuss that relativistic electrons with MeV energies 

are quickly generated in and around the X-type neutral region by utilizing the surfing 

■ acceleration. 

in 
o 



Introduction 



2 ■ Magnetic reconnection is believed to play an important role on plasma heating and acceler- 

■ ation over a broad range of plasma phenomena. Observations of magnetic storms/substroms in 

terrestrial magnetosphere and of solar and stellar flares have shown that enormous hot plasmas 
with temperature of 10^ to 10^ K are generated explosively through magnetic reconnection. In 



these widespread observations, not only the hot thermal plasma population but also nonthermal 
\ particles of MeV energy have been reported to be quickly generated in association with the hot 

plasmas [e.g.. Baker and Stone, 1976; Lin et al, 2003], and the production of energetic particles 
has been recognized as one of the important aspects of reconnection. So far many theoretical and 
observational studies on plasmas heating during magnetic reconnection have been conducted, and 
it is accepted as being beyond doubt that reconnection is one of essential key processes of magnetic 
energy conversion. However, the understanding of the non-thermal particle acceleration mechanism 
remains poor. 

In order to understand the mechanism of particle acceleration, many people have studied the 
charged particle motion under the reconnection electric and magnetic fields, and have conducted 
test particle simulation of the charged particle motion for the prescribed reconnection fields ob- 
tained from the non-linear MHD simulations [e.g., Sato et al, 1982; Scholer and Jamitzky, 1987; 
Birn and Hesse, 1994; Kliem, 1994]. It is discussed that the high energy particle can be effectively 
generated under the reconnection geometry. The acceleration is mainly provided by the direct 
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acceleration of particle with the vx B reconnection electric field. It is also suggested that the mag- 

netohydrodynamic turbulence can contribute to the particle energization as well [e.g., Ambrosiano 
et al, 1988]. Recently Hoshino et al. [2001a; for convenience we shall call this paper as Paper 
1] studied the suprathermal electron acceleration in a particle-in-cell (PIC) simulation where both 
electromagnetic fields and motion of charged particles are self-consistently solved. They found that 
electron acceleration first occurs near the X-type region due to the meandering/Speiser motion 
[Speiser, 1965], and those pre-accelerated electrons arc ftirther accelerated by V B and/or cur- 
vature B drift in the magnetic field pile-up region produced by compression of the reconnection 
outflow. They demonstrated that the spurathcrmal electrons of several 10 keV are generated during 
reconnection in magnetotail, but the acceleration efficiency was not necessarily strong for producing 
MeV electrons. 

Non-adiabatic motion of particle is of major importance in energetic particle acceleration. The 
particle scattering by turbulent electromagnetic waves is recognized as one of the essential processes 
of breaking down the adiabatic motion of particle, and it is the basis of the efficient energy transfer 
process between the particle and the waves. In fact, this kind of scattering process is confirmed 
in the numerical simulation in Paper 1, but it was not a main player of particle acceleration in 
their two-dimensional PIC simulation. Another class of non-adiabatic process can be realized when 
the gyro-radius of the particle becomes comparable to a characteristic scale length of the magnetic 
field curvature [Buchner and Zelenyi, 1989; Delcourt et al, 1996]. In this regime, the particle 
is known to show a chaotic behavior, and the high energy particles can be efficiently accelerated 
by the interaction of the complex particle trajectory under the reconnection electric field. In the 
numerical simulation, this 2nd class of non-adiabatic motion was observed in the pile-up region of 
the reconnecting magnetic field lines where the reconnection outfiow collides with the pre-existing 
plasma sheet plasma. 

It is important to mention that many thin current sheets are often observed during the plasma 
sheet crossing of satellites in the terrestrial magnetotail [Sergeev et al., 1993; Pulkkinen et al., 
1994; Mukai et al., 1998], and the thickness of the current sheet is found to be of order of ion 
inertia length [Fairfield, 1984; McComas et al., 1986; Mitcheel et al, 1990; Asano et al, 2003]. 
Moreover, during the reconnection events with fast plasma flows, a quadrupole magnetic fleld 
structure generated by the Hall electric current and non-gyrotropic ions distributions have been 
reported in the magnetotail observation [e.g., Hoshino et al, 1998; Nagai et al, 2001]. These 
observations suggest that the thickness of the current sheet becomes as thin as the ion inertia 
length just before and during reconnection in magnetotail. 

In such a thin current sheet, energy conversion through kinetic plasma processes plays an 
important role on plasma sheet dynamics and evolution. The micro-scale reconnection structure 
has made remarkable advance lately in both theory /simulation and observational studies [e.g., 
Nishida, 2000; Bij'n et al., 2001]. Thanks to the modern supercomputers, it is now possible to 
discuss the electron dynamics embedded in a macro-scale reconnection structure [e.g., Hesse et al., 
2001; Shinoham et al, 2001; Lapenta et al, 2002; Scholer et al, 2004]. Hoshino et al [2001b] 



-3- 



demonstrated that the strong Langmuiur waves are generated in association with a bump-on tail 
electron distribution function in the plasma sheet boundary. Drake et al. [2003] performed the 
magnetic reconnection simulation with a guide magnetic field parallel to the electric current in the 
plasma sheet, and they observed the intense bipolar electric fields and electron heating near the 
magnetic X-line and separatrices through an electron beam instability. 

Although those intense high-frequency waves observed in the large-scale reconnection simu- 
lations may contribute to the particle acceleration through the wave scattering and/or the direct 
acceleration along the magnetic field, we would like to propose another more efficient acceleration 
process in this paper. We focus on non-thermal electron acceleration by the ambipolar electric fields 
generated around the plasma sheet boundary. Firstly we demonstrate by using the PIC simulation 
that intense ambipolar electric fields are generated under a driven reconnection, namely, in the 
presence of an external driving plasma flow [e.g., Horiuchi and Sato, 1999]. Then we discus that 
the localized ambipolar electric field is particularly important for supra-thermal electron acceler- 
ation through the so-called surfing acceleration [e.g., Sagdeev and Shapiro, 1973; Katsouleas and 
Dawson, 1983]. We argue that relativistic electrons with MeV energies can be quickly generated 
by virtue of the surfing acceleration. 



2. Strong polarization electric fields by driven reconnection 

We study the nonlinear evolution of the collisionless driven magnetic reconnection by using the 
two-dimensional, particle-in-ccU simulation code. The simulation parameters arc basically same as 
those discussed in Paper 1, but we continuously drive a magnetized plasma fiow from both the top 
and bottom lobe regions so as to initiate a fast reconnection [e.g., Horiuchi and Sato, 1999; Bim 
et al., 2005]. The physical system size is L^/X = 38.4 and L^/X = 32, where A is the thickness 
of the initial plasma sheet. We assumed the double periodic boundary in the x and z directions, 
but we set the plasma injection regions at L^/X = ±12, where the magnetized plasma consisting of 
electrons and ions is injected. The external driven electric field at the injection boundary arc given 
by E^ri/^o = ai/cosh^(x/L<^j.i) + 0^2 where L^ri/X = 2 and Eq = VaBq/c. Two cases are simulated: 
one is the weakly driven case with ai = = 0.1, and the other is the strongly driven case with 
= "2 = 0.2. The reconnection structures in the nonlinear saturation stages are almost same, 
but for the weakly driven case it takes twice as long as the nonlinear evolution for the strongly 
driven reconnection. In this paper, we discuss the result for the strongly driven case. 

The numerical grid points are Nx x = 1536 x 1280, and the number of total particles is about 
4 X 10*^. The plasma parameters are as follows: Ion inertia length Va/^i is 1.15 A with 46 grids, 
and electron inertial length c/ujpc is 0.14 A with 6 grids. The mass ratio of ions to electrons mi/rrif, 
is 64, for simplicity. The ratio of the speed of light to the Alfven velocity c/Va is 16.8. The plasma 
temperatures are set as follows: Ti^pg/Te^ps = 4, Ti^obe/^i.ps = 1/25, and Ti^iobe = ^e^iobe, where 
the subscript "i" means ion, "e" electron, "ps" the plasma sheet, and "lobe" the lobe, respectively. 
The ratio of the Alfven velocity Va to the ion thermal velocity vth,i = \/2I|^^^7^ is 1.15, where 
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the Alfven velocity is defined by the lobe magnetic field Bq and the plasma sheet density uq. The 
initial lobe plasma density is set to be 5 % of the plasma sheet density no- 

Figure 1 shows a snapshot of the nonlinear stage of a driven reconnection at t/r^ = 52.2, where 
TA = ^/Va is the Alfven transit time. In the top panel, the color contour of the electric field is 
shown together with the magnetic field lines denoted by the solid lines, while the bottom panel shows 
the electric field E^. We find that a pair of V-shaped strong electric fields are induced almost along 
the plasma sheet boundary, namely around the boundary between the lobe and the plasma sheet. 
The electric fields point towards the plasma sheet, which directions are almost perpendicular to the 
local magnetic field line. These ambipolar electric fields are known to be produced in association 
with Hall electric current in a thin current sheet due to the inertia difference between ions and 
electrons, when the thickness of the plasma sheet becomes comparable or less than the ion inertia 
length [Hoh, 1966; Hoshino, 1987]. For the thin plasma sheet, ions are unmagnetized and do not 
move with the magnetic field line motion, while electrons remain magnetized and are frozen in the 
magnetic field line. Therefore, the electron density becomes higher than the ion density in the 
thin plasma sheet, and an ambipolar /polarization electric field is produced in order to maintain 
the charge neutral plasma. In non-linear stage of reconnection, those electric fields are localized in 
the boundary region. From our simulation results, we find that the amplitude of the polarization 
electric fields reaches up to about 10 - 20 times of the driven electric field E^i^i, and the typical 
region width where the polarization electric fields are induced is less than the ion inertia scale. 

In addition to the above large scale polarization electric field, we also find a series of large 
amplitude electrostatic wave emission propagating outward from the X-type region in the boundary 
region. In the spontaneous reconnection in our previous study [Hoshino et ai, 2001b], we found 
a signature of the Langmuir wave emission, but the amplitude remained too small to recognize it 
as a series of wave trains. The large amplitude waves in Figure 1 can be produced as the result 
of the nonlinear evolution of the Langmuir waves in the driven reconnection. It is interesting to 
note that Drake et al. [2003] reported the similar bi-polar electric field generation near the X-type 
region with a guide magnetic field parallel to the electric current, namely, the situation for the 
magnetopause reconnection. 

In the bottom panel of Figure 1, we find that E^ is enhanced for both top and bottom regions 
around {X,Z) = (0, ±4). The top region has the positive Ez, while the bottom is the negative 
one. These regions with a finite electric field Ez correspond to the front side of the compressed 
plasma propagating from the driven boundary at Z/\ = ±12. We would like to mention that the 
time evolution of reconnection was speeded up when the location of the driven boundary is set 
to be closer to the neutral sheet, but in order to get the better separation between the X-type 
reconnection structure and the above front side region of the driven flow, we chose the boundary 
location of the external magnetized plasma flow at Lz/\ = ±12. 

It is also noteworthy that the external driven system dramatically changes the time evolution 
of system [e.g., Hayashi and Sato, 1978; Hoshino, 1991; Pritchett and Coroniti, 1995; Horiuchi and 
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Sato, 1999; Birn et al, 2005]. In the driven current sheet, the current layer becomes thin with 
increasing the lobe magnetic field. The time evolution of the electric current at the neutral sheet 
Jns for one-dimensional slab geometry can be expressed by [Drake et al., 1981; Hoshino, 2004] 

„.(0)(^)*'\ (1) 

where B^ri is the amplitude of magnetic field at the driven boundary, and F is the ratio of specific 
heat. Then the time evolution of thickness of the driven current sheet Xdn on B^ri/Jus can be 
evaluated as, 

Ao ^Wa^)) ' 
where Aq is the initial current sheet thickness. For the case of one-dimensional plasma sheet 
thinning, F may be taken to be 2 because of two-dimensional velocity space compression. However, 
if the time scale of the isotropization of the anisotropic velocity distribution function is faster than 
the plasma compression time scale, we can use F = 5/3. If we could assume an isothermal state 
with F = 1, X(iri oc the thickness of the current sheet is crucially affected by the boundary 

magnetic field. The reason why the strong influence appears for an isothermal state is that the 
thermal energy is extracted due to the heat flux escape from the neutral sheet, and that in order 
to maintain the pressure balance the current sheet quickly shrinks. For the case of non-uniform 
external driven boundary condition such as E^ri oc cosh~^(x), it is easily expected that the heated 
plasma at the center of the highest density region can be removed along the magnetic field line. 
Having those effects in our mind, we may expect that the driven force in the reconnection setting 
has a strong influence on the current sheet thinning even before the onset of reconnection. 



3. Nonthermal electrons and surfing acceleration 

Shown in Figure 2 is energy spectra of electrons integrated over pitch angle in the whole 
reconnection region. The energy e is defined as e = (7 — l)meC^, where 7 is the Lorentz factor of 
the particle. Two dashed curves show the spectra at an early stage of reconnection at I/ta = 7.5 
and the nonlinear stage at t/r^ = 44.7 before the onset of the strong V-shaped polarization electric 
field seen in Figure 1. The spectra with the solid curves are about the well-developed nonlinear 
stages at I/ta = 52.2 and I/ta = 59.6 when the electrostatic wave emission becomes strongly active. 
As the time goes on, one can find the relativistic electrons are generated, and a large amount of 
relativistic electrons are quickly produced, and the nonthermal energy density normalized by the 
thermal energy one becomes 10%. 

In the spontaneous reconnection discussed in Paper 1, we already found the production of 
suprathermal electrons, but the nonthermal electron energy density was about 1% compared with 
the thermal energy density. We will discuss later that the difference of energization is due to 
the amplitude of a pair of the V-shaped polarization electric fields excited around the boundary 
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between the lobe and the plasma sheet. In fact, the amplitude of the polarization electric field in 
the spontaneous reconnection is about E/Eq ~ 1, while that in the driven reconnection is about 
E/Eo - 5 - 10. 

In order to understand how and where the high energy electrons are generated, we investigate 
an electron trajectory. Shown in Figure 3 is a typical trajectory of electron that is strongly accel- 
erated during the interaction with the plasma sheet. The solid curve is the trajectory in (X, Z) 
space, and the dashed curve shows the energy history as a function of X. The trajectory starts at 
(X, Z) ~ (8.7, 2.3) denoted by the position 1, where the electron is situated in the lobe region, and 
the electron remains a very low energy. As the electron is moving toward the X-type neutral line, 
it gradually gains the energy of e/rrieC? ^ 0.3 at (X, Z) ~ (1, 0) indicated by the position 2. During 
the Speiser motion around the X-type neutral line, the electron moves from the position 2 to 3, 
and the electron gains further energy up to e/mgC^ ~ 0.6. After the electron is ejected outward 
from the X-type neutral region, it travels toward the magnetic field line pile-up region, which may 
correspond to the region between the position 4 and 5. In this pile- up region, V B/curvature B drift 
motion plays an important role on particle energization, and we find that the electron is accelerated 
up to e/mgC^ ~ 0.7 during a chaotic motion around the position 5. 

We discussed two-step acceleration mechanism of electron in Paper 1: one is the standard 
Speiser /meandering motion around the X-type region, and the other is V B/curvature B drift 
acceleration for non-adiabatic particles with k ~ 1 in the magnetic field line pile-up region, where 
is the ratio of the curvature for the magnetic field and the Larmor radius [Buchner and Zelenyi, 
1989]. In the driven reconnection system, the reconnection fields grow up much faster than the 
spontaneous reconnection, and the point of this paper is that we found the pre-stage acceleration 
before the Speiser acceleration. The acceleration seems to be switched-on during the nonlinear 
evolution of the strong electric fields shown in Figure 1. 

For convenience we shall call the pre-acceleration as the 1st stage acceleration, the Speiser 
acceleration as the 2nd stage, and the acceleration in the pile-up magnetic field region as the 3rd 
stage acceleration. We think that all three acceleration phases are important for the acceleration 
of electrons to high energies. Figure 4 summarizes the above acceleration processes. Except for 
the pre-acceleration of electron during the cold electron transport stage from the lobe toward the 
X-type neutral line, the acceleration processes listed here are same as those discussed in Paper 1. 

Let us discuss the 1st stage acceleration in detail. We think that the acceleration mechanism 
in the 1st stage is understood by a "surfing acceleration" due to the polarization electric field 
induced in the boundary between the lobe and the plasma sheet, because the acceleration occurs 
in the direction perpendicular to both the magnetic field and the polarization electric field. Let us 
explain what the surfing acceleration is and how the surfing acceleration works in the reconnection 
region. Shown in Figure 5 is a cartoon of the key elements of the electromagnetic fields in the surfing 
acceleration. The dashed lines show the magnetic field lines around the reconnection region, and we 
take the tilted coordinate system of {X', Z') where the X' axis is parallel to the local magnetic field. 
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The polarization electric fields points the negative Z' direction almost perpendicular to the local 
magnetic field. The equation of motion for an electron in this coordinate system can be written as 

eE^, ~ 0, (3) 

eEy - -Vz'B^, ~ -eEy, (4) 

e 

eE^i ^ -VyB:^,. (5) 

Wc have assumed that the field-aligned electric field E^i is negligible for the surfing acceleration, 
but note that the small parallel electric field may plan an important role on the thin current sheet 
structure [Zelenyi et al., 2004]. In Eq.(4) we have neglected the Lorentz force term of v^'B^', 
because this term remains small after averaging it during the trapping motion in the z-direction. 
Let us also assume that the energy of the electron transported from the lobe is smaller than the 
potential energy ecj), where E^i = —d(f)/dz'. Then the electron can be trapped by the potential well, 
and we assume that v^i is very small for the early phase. During this trapping phase, the electron 
can be accelerated by the reconnection electric field Ey, if the electric force of eE^i can balance 
with the Lorentz force of eVyB^' /c. 

If the Lorentz force becomes larger than the electric force during this acceleration, the particle 
can escape from the barrier of the polarization electric field E^', and the energization of particle 
will stop. The maximum velocity of the surfing acceleration can be roughly estimated by 

^^max ~ C— . (6) 
Bx' 

One of the key points of the surfing acceleration is that the energization of particles comes from 
the resonance with the Ey electric field, even though the final velocity Vmax is expressed by E^^/. 
Another point is the localization of the E^' electric field. If E^' were uniform in space, the surfing 
acceleration is replaced by the standard E x B drift. In the E x B drift case, the energization 
during the initial pick-up phase can be provided from the Ez' electric field under Ez' 3> Ey. 

The other interesting behavior of the surfing acceleration appears under a strong Ez' field. 
If the electric field E^' is larger than the magnetic field B^', the particle cannot escape from the 
barrier of the potential well, and it will be able to gain quickly a relativistic energy. This kind of 
acceleration does not appear for most of MHD phenomena where E/B ^ 0{Va/c), but this is not 
the case for the driven reconnection in a thin plasma sheet, because the localized electric field can 
become as large as the reconnection magnetic field near an X-type region. 

Figure 6 is the velocity hodogram in (Py,Pz) for the same electron shown in Figure 3 (where 
Pj = jVj/c), and the same time interval is chosen. The markers of positions 1 to 5 are chosen 
to be the same time stage as Figure 3 as well. In the initial phase before 1st stage acceleration, 
the electron is located around {Py,Pz) = (0.0,0.0). During the 1st stage acceleration, the electron 
moves toward position 2 with {Py,Pz) = (0.8,-0.2) along the Py velocity axis, which suggests 



dpx' 

dt 

dpy 

dt 
dpz' 

dt 
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that the particle gains energy by traveling toward the anti-direction against the electric field Ey. 
This behavior confirms that the DC acceleration occurs in a finite magnetic field region, which 
is known as the surfing acceleration [e.g., Katsouleas and Dawson, 1983; Hoshino and Shimada, 
2002]. In 2nd stage of meandering/Speiser acceleration, the electron further moves to position 3 
with {Py,Pz) = (1.0,-0.2), and after that the electron changes the velocity direction due to the 
Lorentz force and is ejected from the X-point region. After position 5, the large circular motion 
shows the gyro-motion of the accelerated electron. 

If the gyro-radius for the energetic particles generated during 1st and 2nd acceleration stages 
becomes comparable to the local curvature of the magnetic field pile-up region, the electron motion 
becomes non-adiabatic so that some of electrons can be further accelerated under the reconnecting 
magnetic fields with the mirror geometry. Prom our PIC simulations, we conclude that together 
with the meandering/Speiser acceleration, the surfing acceleration can contribute to the production 
of energetic electrons whose gyro-radii are comparable to a local magnetic field curvature. 



4. Conclusions and Discussion 

We demonstrated for the first time the electron surfing acceleration during reconnection in 
a thin plasma sheet by using the PIC simulation. However, we used a reduced mass ratio of 
m,i/me = 64 in order to save a computational CPU time, and it is necessary to discuss carefully 
the validity of the surfing acceleration in more general plasma parameters. 

One of the main points in the surfing acceleration is the amplitude of the polarization electric 
field in association with the Hall electric current, because the efficiency of the surfing acceleration 
is controlled by the amplitude of V-shaped polarization electric field induced around the X-point 
region. Let us quickly estimate the amplitude of the electric field. The electrostatic potential 
energy induced in the boundary may be equated to the ion thermal energy. The ions are assumed 
to be in a Boltzmann-Maxwell state, because ions are unmagnetized in a thin plasma sheet, and 
we roughly obtain, 

e<j) ~ Ti. (7) 

The thickness I of the electrostatic potential may be estimated as an intermediate scale between 
the ion inertia and the electron inertia length from our simulation result, and for a moment we may 
assume that I is estimated by the hybrid inertia scale, i.e., I = c/ y^UJ^^oJ^. Then the polarization 
electric field Ep can be obtained as, 

^ ^ A ^ V ^ ^ ( 

Bo IBq \UJpeJ \l.^UJpeOJpiJ \mi 

where ujce = eBo/nieC is defined by using the lobe magnetic field, and Upe is for the plasma sheet. 
We have assumed that the ion thermal velocity is almost equal to the Alfven speed Va- The ratio 
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of ujce/^pe is a Small value in the magnetotail. For a typical tail parameters, 



u)rf. „ „ . o f B(\ \ f 0.5cm 



V20nT 



3\ 1/2 



^~8.8xlO-M^ pi:^ , (9) 



no 



and we will get the amplitude of the polarization electric field in magnetotail is 

- -o(^Vf^V'7^Vv/. (10, 

~ 13 
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^0.5cm ' 
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We find that the polarization electric field ranges from 10 mV/m (for I = c/upi) to 500 mV/m (for 
I = c/upe). 

If we normalize the porlarized electric field Ep by the typical MHD electric field of BqVa/c, 
we get 

(12) 



BqVa/c ^me\ujpelj \uJpil , 

Even if the scale length I is equal to the ion inertia length of c/u)pi, the polarization electric field 
Ep becomes of order of the typical MHD electric field of Eq = BqVa/c, which is still larger than 
the reconnection electric field E^ = arEo in the order of 0.1 x BqVa/c (~ 1 mV/m in magnetotail), 
where we have assumed the reconnection rate ~ 0.1. 

Next let us estimate the maximum electron energy Smax during the surfing acceleration. Prom 
Eq.(4), we obatin, 

£max ~ / eE,Vy{t)dt ~ V^^^* ' (13) 

where T* is the time interval of the surfing phase. We have assumed that the initial electron energy 
is negligible compared to the final energy £max- The time inverval of the surfing phase T* may be 
estimated as, 

~ -— , (14) 

ar VeO 

where L is the thickness of the plasma sheet where the polarization electric field is induced, and 
the distance along the magnetic field from the lobe to the X-type neutral point can be multiplied 
by the geometrical fact of l/a^ . VeO is a typical electron velocity along the plasma sheet boundary 
during the surfing phase, and it may have an intermediate velocity between the lobe and plasma 
sheet thermal velocities. Then we obtain, 



rUeC^ 2 \UJpeJ \ C J \ VeO 



0(10-2), (15) 



where Tj/Tg = 9, which is known to be about 5 in the terrestrial magnetotail. We have assumed 
that Vth,e/veO is of order of 0.1 ~ 1. The above equation suggests that electrons can be accelerated 
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up to several keV during the surfing acceleration, therefore those electrons can satisfy non-adiabatic 
condition with k ~ 1 for the reconnecting magnetic field line in a thin plasma sheet [Hoshino et 
al, 2001a]. Thanks to this pre-acceleration during the plasma transport from the lobe toward 
the X-type region, V B and/or curvature drift as the 3rd stage of acceleration can give further 
energization (see Figure 4). We think that, by combining all three acceleration steps, electrons 
with MeV energies can be quickly produced during magnetic reconnection. 

It should be noted that the estimation of the maximum energy in Eq.(15) is valid only when 
the trapping condition of electrons with E^' ~ VyB^' /c is satisfied during the surfing acceleration at 
the plasma sheet boundary. As the electron travels towards the X-type region, it is expected that 
the amplitude of the local magnetic field Bj.i will decrease. If so, the trapping condition could be 
still satisfied for a large Vy, and the electron can gain energy until the maximum energy estimated 
by Eq.(15). More careful analysis of the trapping condition would be necessary in order to clarify 
the efficiency of the surfing acceleration. 

It might be interesting to note the structure of the non-thermal energy spectrum. Shown in 
Figure 7 is the energy spectrum at I/ta = 59.6, which is the same energy spectrum in Figure 
2 except for the double logarithmic scales. We find that the spectrum can be approximated by 
a double power-law spectrum, in which the lower energy band has the spectral index of about 3 
and the higher band with about 7. The brake energy of the double power-law spectrum is about 
^max/fneC^ — 0.7. In our simulation parameters, we used 9 = A, (ujce/^pe)'^ = 0.05, and Lujpi/c ~ 2, 
and the typical electron velocity v^o ~ ^.f)Vth,e- Therefore, from Eq.(15), the maximum energy 
expected from the surfing acceleration is estimated as £rnax / '^e(? — 0.8. Since the maximum energy 
of the surfing acceleration is very close to the simulation result, we think that the lower-energy 
population of the double power-law energy spectrum can be generated by the surfing acceleration, 
and the higher-energy population is mainly produced by V B/curvature B drift motion. It would 
be interesting to mention that the double power-law spectrum is observed in the solar corona by 
the RHESSI satellite [e.g., Lin et al., 2003], and the kappa distribution functions are often observed 
in the terrestrial magnetotail [e.g., Christon et al., 1989]. The surfing acceleration might be one of 
mechanisms to explain the observed power-law spectrum. 

Finally, we should like to make the point that it may be possible to explore many aspects of the 
electron surfing acceleration from in situ observations in the terrestrial magnetosphere. One of the 
key points of the surfing reconnection is to confirm observationally a structure of polarization electric 
field. There are several observations regarding to the spiky electric field structure in the boundary 
between the lobe and the plasma sheet [e.g., Cattell et al., 1994], but so far those observations 
are not discussed in the framework of the model presented here. Further investigations on the 
small-scale electric field structure could give us better understanding on the electron acceleration 
in the course of magnetic reconnection. 
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Fig. 1. — The magnetic field lines and the color contours of the electric fields at I/ta = 52.5 are 
shown. Both Ex and Ez are normalized by ii^o = ^ Bq/c. Large amplitude polarization electric 
fields induced around the boundary between the lobe and plasma sheet. A series of bipolar electric 
field structure propagating outward from the X-type region along the magnetic field lines can be 
also seen. 
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Fig. 2. — Electron energy spectra integrated over pitch angle in the simulation domain at I/ta = 
7.5,44.7,52.2, and 59.6. 
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Fig. 3. — A typical electron trajectory with the surfing acceleration in the x — z plane (solid line), 
and its energy history as a function of x (dashed line). 
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X 

Fig. 4. — Sketch of acceleration processes during magnetic reconnection. Ep is the polarization 
electric field induced in a thin plasma sheet, which plays an important role on surfing acceleration. 
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Fig. 5. — Illustration of the electric and magnetic field structure for electron surfing acceleration 
in the reconnection region. 




Fig. 6.— 



A velocity hodogram for the same electron analyzed in Figure 3. 




Fig. 7. — Double power-law energy spectra at t/TA = 59.6 in the double-logarithmic scale. 



